Form  Approved 
OUB  Ho.  OT04-01M 


AD-A282  701 


>N  PAGE 


1.  AGENCY  USE  ONLY  (Loavo  blank) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  UA  l  e 

March  1987 


PEGS  Long  Conduction  Time  Theory 


1  Imw  p«r  rasponss.  Including  Ilia  Una  tai  reviewing  instructions,  searching  existing  data  sources. 
Ueetion  al  Information,  tend  comments  regarding  this  burden  estimate  er  any  ether  aspect  et  this 
Ungten  Headquarters  Services,  Directorate  ter  Inferatatten  Operations  and  Itaperts,  if  If  Jo  Her  sen 
ament  and  Budget,  Paperwork  Reduction  Project  (0704-01  ft).  Washington.  DC  tOSOS. 


t I  3.  REPORT  TYPE  AND  DATES  COVERED 

I  Technical  - 


S.  FUNDING  NUMBERS 

DNA-001-85-C-005 1 


6.  AUTHOR(S) 

E.  Waisman,  D.  Parks,  P.  Steen,  A.  Wilson 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

S-Cubed 

A  Division  of  Maxwell  Laboratories,  Inc. 

P.O.  Box  1620 
La  Jolla,  CA  92038-1620 


S.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


SSS-DFR-89-10651 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 


This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under 
RDT&E  RMSS  Code  B3230  85466  G37TMXGE  00018  25904  D. 


12a.  DISTRIBUTION/ AVAIL  ABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Public  Availability 


13.  ABSTRACT  (Maximum  200  words) 


PEOS  technology  scaling  to  the  long  conduction  time  (~>  200ns)  and  10-20  MA  current 
levels  required  by  the  10X  machine  is  discussed.  The  approach  taken  is  to  employ  the 
results  of  the  electromagnetic  numerical  modeling  we  have  performed  for  conduction 
times  of  50  to  80  ns  in  constructing  the  functional  dependence  of  the  switch  operation 
on  its  physical  parameters.  BLACKJACK  5  and  ACE  PEOS  experiments  as  well  as  results 
from  fully  electromagnetic  particle  in  cell  simulations  are  compared  with  the  heuristic 
ion  piston  model.  Implications  of  this  analysis  on  switch  length,  radius  and  density  for 
the  10X  machine  are  given. 


15.  NUMBER  OF  PAGES 
21 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  I  IS.  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION!  20.  LIMITATION  OF  ABSTRACT 


OF  REPORT 

UNCLASSIFIED 


NSN  7540-01  -280-5500 


OF  THIS  PAGE 

UNCLASSIFIED 


OF  ABSTRACT 

UNCLASSIFIED 


brio  i 


Standard  Form  290  (Rev.  2-89) 
Preserved  by  ANSI  Std.  Z39-18 
298-102 


SSS-DFR-89-10651 
W/O  11283 


PEOS  LONG  CONDUCTION  TIME  THEORY 


E.  Waisman 
D.  Parks 
P.  Steen 
A.  Wilson 


S-CUBED 

A  Division  of  f  axwell  Laboratories,  Inc. 
P.O.  Box  1620 
San  Diego,  CA  92038-1620 


Accesion  For 

NTIS  CRA&I 
DTIC  TAB 
Unannounced  □ 
Justification 


By . 

Distribution  / 


Availability  Codes 


Dist 

Avail  and/or 

Special 

/?'/ 

Contract  No.  DNA-001-85-C-0051 

This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under 
RDT&E  RMSS  Code  B3230  85466  G37TMXGE  00018  25904  D 


Presented  at: 

DNA  Pulsed  Power  Review  and  SDIO  Update 
Las  Vegas,  Nevada 


94-23362 


31  March  -  2  April  1987 


94  7  25  0  07 


PEOS  LONG  CONDUCTION  TIME  THEORY 


E.  Waisman,  D.  Parks,  P.  Steen,  and  A.  Wilson 
VG  I  General  program  overview. 

VG  2  The  objective  is  to  discuss  how  PEOS  technology  scales  up  to  the  long 
conduction  times  (j>  200  ns),  and  10  to  20  MA  currents  required  by  the  10X 
machine. 

The  approach  is  to  employ  the  results  of  the  electromagnetic  numerical 
modeling  we  have  performed  for  conduction  times  of  50  to  80  ns,  in 
constructing  the  functional  dependence  of  the  switch  operation  on  its  physical 
parameters. 

We  compare  the  BLACKJACK  5  and  ACE  PEOS  experiments,  as  well  as  the 
E&M  simulations,  with  a  heuristic  model:  the  ion  piston  mechanism. 

Implications  of  this  analysis  on  switch  length,  radius  and  density  for  the  10X 
machine  are  given. 

VG  3  A  schematic  representation  of  the  particle  electromagnetic  model  is  shown. 

Ions  and  electron  self-consistent  trajectories  are  followed  in  a  planar  mesh  of 
length  >  £,  and  height  DR,  where  £  is  the  part  of  the  computational  grid 
occupied  at  t=0  by  the  initial  carbon  neutral  plasma,  of  ion  density  N  and 
injection  velocity  v.  DR  represents  the  A-K  gap.  The  total  length  of  the 
mesh  includes,  besides  £,  sections  to  the  left  and  right  of  it  initially  empty  of 
particles. 

The  input  (upstream)  side  is  driven  by  Vqc,  the  open  circuit  voltage  of 
BLACKJACK  5  shot  #1057  in  most  cases,  coupled  by  a  L-Z  circuit  with  L 
and  Z  constant  in  time  (upstream  inductance  and  line  impedance,  respectively). 
The  output  side  (downstream)  is  a  L-Z  circuit,  where  L  is  the  downstream 
inductance  to  the  load,  and  Z  is  the  time  dependent  magnetic  limited 
impedance  of  a  ring  e-beam  diode. 
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VG  4 


Shown  in  this  viewgraph  is  the  comparison  with  the  experiment  of  the 
electromagnetic  simulation  using  the  actual  length  and  gap  of  shot  PEOS 
BLACKJACK  5  1057,  i.e.,  Z  =  18.1  cm  and  DR  =  4.2  cm.  The  density  and 
injection  velocity  of  the  C  plasma  in  the  simulation  were  N  =  4  x.  10  /cm 
and  Vp  =  10  cm//is,  respectively.  This  density,  for  the  given  switch  size,  is 
near  the  limits  of  what  these  calculations  can  treat  in  reasonable  computing 
times,  even  at  relatively  low  resolution. 

Agreement  between  the  E&M  simulation  and  experiment  is  semiquantitative 
(we  remark  that  the  actual  density  distribution  is  only  known  to  order  of 
magnitude  accuracies  in  the  experiments).  The  simulation  shows,  for  this 
density  of  4  x  10  /cm  of  C  ,  shorter  conduction  time  than  the  observed 
one,  and  the  presence  of  a  "numerical”  foot  for  the  load  current.  No  evidence 
of  loss  of  magnetic  insulation  or  sizable  switch  current  after  opening  is 
observed  in  the  simulation.  Perhaps  this  difference  between  simulation  and 
experiment  is  due  to:  (i)  In  the  calculations  no  geometric  details  such  as 
corners,  convolute  posts,  etc.  in  the  vacuum  section  between  switch  and  local 
are  considered,  (ii)  We  used  a  simplified  impedance  model  which,  upon  the 
load  getting  to  a  voltage  ~  100  kV,  goes  from  8  0  to  about  less  than  1  fl  in 
5  ns.  The  experimental  observation  shows  Z  ~  1.5  to  2  fl  at  that  stage. 

At  this  point  of  the  talk  a  video  cassette  was  shown  of  ion  density,  electron 
density,  and  magnetic  field  profiles,  obtained  from  the  E&M  calculation  for  a 
switch  of  length  £  =  9  cm  and  with  all  other  parameters  the  same  as  in  the 
previous  viewgraph.  The  images  show  the  formation  of  a  current  channel  at 
early  times  and  the  formation  of  an  ion  piston  (snowplow)  at  times  before 
opening. 

VG  5  Shown  in  this  viewgraph  is  the  definition  of  T£  and  Tq,  conduction  and 
opening  time,  respectively,  for  the  E&M  simulations. 

Such  definition  is  necessary  because  we  observe  a  “numerical”  foot  in  our 
calculated  load  currents,  that  we  conjecture  is  due  to  spurious  numerical 
diffusion,  consequence  of  discretization  errors. 
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VG  6 


Numerical  scaling  studies  were  done  by  keeping  Vqc,  total  inductance 

(upstream  +  switch  +  downstream)  constant,  and  using  the  same  load  model. 

The  plasma  length,  density,  and  gap  were  changed,  as  indicated  in  the 

viewgraph;  also  R  (R^q^  +  ^c^thode^  was  vari*d  **  since  our  geometry  is 

planar,  R  only  enters  in  the  relationship  between  current  and  magnetic  field  in 

the  circuit-boundary  conditions  Also  the  injection  velocity  was  changed  from 
6  7 

10  to  10  cm/s,  and  we  consider  two  cases,  identical  in  ail  other  respects, 
differing  by  charge  state,  i.e.,  C~*~  and  C In  total  about  20  different  cases 
were  run  for  times  between  30  and  80  ns. 

VG  7  Scaling  of  the  numerical  simulations  follow  the  indicated  functional  dependence. 
Later  in  the  sequence  of  viewgraphs  we  explain  the  ion  piston  model. 

VG  8  This  viewgraph  shows  schematically  the  differences  in  load  turn-on  between 
theory  and  experiment. 

VG  9  In  this  viewgraph  the  formulation  of  numerical  diffusion  is  introduced.  It  is 

shown  that  the  “numerical  diffusion  time”  is  proportional  to  the  number  of 

9  10 

macroparticles  (typically  10  to  10  electrons  or  ions)  per  computational  cell. 
These  effects  seem  to  be  bothersome,  but  not  to  overwhelm  the  physics 
studied:  for  the  case  shown  before,  simulation  of  shot  #1057,  the  diffusion 
time  is  of  the  order  indicated  in  this  viewgraph,  i.e.,  somewhat  larger  than  the 
calculated  conduction  time. 

VG  10  The  equations  describing  a  1-D  leaky  ion  piston  are  introduced.  Here  p  is 
mass/unit  of  the  plasma  area,  p  its  volume  mass  density,  x  is  the  position  of 
the  ion  piston  front;  x  =  0  at  the  initial  edge  of  the  switch,  M  is  the  ion 
mass.  The  “leakiness”  is  represented  in  equation  (2)  in  the  term  p/T,  r  is 
the  characteristic  time  that  an  ion  of  mass  M  would  take  to  traverse  the 
colissionless  length  c/w^,  if  accelerated  by  the  force  eBx/c,  where  x  is  the 
velocity  of  the  ion  piston  front.  Equation  (4)  is  the  circuit  equation  valid  for 
the  conduction  phase,  i.e.,  the  time  it  takes  for  the  ion  piston  to  reach  the 
end  of  the  switch,  that  is  to  go  from  x=0  to  x=2;  g  is  the  A-K  gap. 

VG  11  The  predictions  of  the  ion  piston  model  are  shown  in  this  viewgraph.  With 

T  the  conduction  time,  B(T  )  the  magnetic  field  at  that  time  (upstream,  since 
c  c 

for  t  S  Tc  the  load  current  is  zero),  and  2  the  switch  length. 
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The  observations  of  ACE  and  BLACKJACK  5  in  comparable  geometry  (same 
cathode  radius,  same  length  and  same  plasma  timing)  seem  to  obey  this 
prediction.  We  remark  that  this  prediction  about  BTc  agrees  much  better 
with  the  observations  than  the  threshold  current  model,  which  would  have 
predicted  equal  currents  for  both  cases,  irrespective  of  the  factor  of  three  in 
conduction-time  of  ACE  with  respect  to  BLACKJACK  5. 

VG  12  The  NRL  model  scaling  of  the  peak  switch  magnetic  field  (a  I  ft  )  with 

sw  c 

switch  length  L  is  compared  with  BLACKJACK  5  and  the  recent  ACE 

8W 

PEOS  experimental  data.  The  linear  scaling  predicted  by  that  model  is  shown 
to  be  poor. 

VG  13  The  scaling  of  magnetic  Held  (I  / r)  with  the  average  ion  piston  velocity 
(L  /r)  where  T  is  the  conduction  time  is  shown.  The  prediction  of  the 
model  is  good,  indicating  that  the  model  can  be  utilized  for  extrapolating  to 
the  long  conduction  times  needed  in  the  10X  machine. 

VG  14  This  viewgraph  illustrates  a  comparison  of  the  upstream  current  vs  time 

between  experiment,  shot  #1235,  the  ion  piston  model  and  a  constant 

inductance  assumption  for  the  duration  of  the  conduction  phase.  The 

experimental  load  current  is  shown  to  indicate  when  the  conduction  phase  is 

12  3 

over  -  turn-on  of  load.  The  density  shown  Ncar^oQ  =  1.04  x  10  /cm  ,  is 
the  one  which  gives  the  best  fit  to  the  experiment  using  the  ion  piston  model, 
such  as  to  have  the  ion  piston  front  reach  the  downstream  end  of  the  plasma 
at  T£.  L-DOT  =  dL/dt  effects  are  seen  to  be  significant. 

VG  15  Conclusions  are  self-explanatory. 
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PLASMA  EROSION  OPENING  SWITCH 


OVERVIEW 
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Viewgraph  2 


SCHEMATICS  OF  LUMP  CIRCUIT  FOR  BLACKJACK  5 
PEOS  ELECTROMAGNETIC  PARTICLE  SIMULATIONS 


BLACKJACK  5  PEOS  SHOT  1057 


Viewgrapi 


DEFINITION  OF  T-  AND  T, 


Viewgraph  5 


NUMERICAL  SCALING  STUDIES 


Viewgraph  6 


NUMERICAL  SCALING  RESULTS 
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LACK  OF  AGREEMENT  WITH  ION  PISTON  MODEL  IN  ABSOLUTE  TERMS 


THEORY  EXPERIMENT  DIFFERENCES 


Viewgraph  8 


NONPHYSICAL  NUMERICAL  EFFECT 


Viewgraph  9 


SCALING  PREDICTED  BY  ION  PISTON  MODEL 
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EXPERIMENT  —  NRL  MODEL  COMPARISON 
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EXPERIMENT  —  ION  PISTON  MODEL  COMPARISON 


17,894  Viewgraph  13 


Viewgraph  14 


CONCLUSIONS 
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FURTHER  VALIDATION  OF  SCALING  WITH  N,  2  AND  SWITCH  RADIUS  IS 
ESSENTIAL 
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